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Abstract — In the mathematical discipline of linear algebra 
there are many operations defined on matrices such as addition , 
multiplication and matrix inverse etc. Actually we note that in 
the real world it is not feasible for most of the matrix 
computations to be calculated in an optimal explicit way, such as 
matrix inversion . Matrix decomposition is a fundamental theme 
in linear algebra and applied statistics which has both scientific 
and engineering significance. Now suppose we are given a 
matrix A and we have to transform it into a product of two 
matrices ,say we want to write A=BC therefore in this case a 
numerically stable and fairly fast scheme is desired to compute 
the two matrices A and B. In general we mean by the sentence 
matrix decomposition or matrix factorization the factorization 
of a matrix into a product of two or more matrices, so the 
process of converting or transforming a a given matrix in one 
side to product of two or more matrices in the other side is 
called matrix decomposition or matrix factorization. Really 
matrix factorization or matrix decomposition has numerous 
algorithm like LU factorization ,QR decomposition ,rank 
factorization, singular value factorization etc . But most of these 
algorithms are complicated ,so our goal here is to find a simple 
algorithm with which we can factorize or decompose the 
specified matrix and this was done by making use of the inverse 
of the matrix . Sometimes we need to factorize a matrix 
according to different constraints like non- negative matrix 
factorization and this happened if we have a non-negative 
matrix and we want to factorize it also into non- negative ones, 
or any other constraints due to the nature of the problem which 
was represented by the matrices ,some methods of 
decomposition factorizes approximately not exactly but here in 
our easy new approach the factorization will be exact . 

In this paper we will introduce a very simple algorithm with 
which we can factorize any matrix exactly to a product of any 
numbers of matrices by factorizing again each matrix of 
matrices which form the product . 

Index Terms — Invertible matrix, matrix decomposition, 
noisy matrix, the entries of the matrix, singular matrix, 
constraints , matrix multiplication , associative law, the rank of 
the matrix. 


I. Introduction 

This paper propose an easy new method for matrix 
decomposition and it contains four sections, in section 2 there 
is a theorem with its proof using the notion of the determinant 
and of the matrix and its properties, in section 3the invertible 
matrix play an important role for matrix decomposition 
besides the associative law of matrix multiplication, it also 
contains some numerical example to make the idea very clear, 
section 4 is the conclusion of the paper. 

We note matrix factorization refers to transformation of a 
given matrix into a given multiplication form [1], 

Now let matrix A be the product of the matrices W and H, 
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A=WH 

Matrix multiplication can be implemented as computing the 
columns vectors of A as linear combinations of the columns 
vector in W using coefficients supplied by the columns of H, 
so that each column of A can be computed as follows 
a, - w l h I t 

Where a j is the ith column vector of the matrix A and /? is 

the ith column vector of the matrix H . 

To make it clear we will give a common example which is 
facial image decomposition^]. 

Let A a matrix with a real-valued, nonnegative entries of 
size m x n such that 
A = WH 

Where W is m x r , H is r x n , and both factors have none 
negative entries . If the columns of A are pixels of a facial 
image, the columns of W may be facial features such as eyes 
or ears , and the coefficients in H represents the intensity of 
these features . Non-negative matrix factorization is used in 
abroad range of applications including graph clustering [3], 
Protein sequence motif discovery [4], and hyper spectra 
unmixing [5] . 

For factor analyses which is a model for real- valued data we 
can construct an equation of four matrices in the form 
A = WH + E w here is E called a noisy matrix or the error 
matrix .and we can use the least square method to find the 
minimum error so as to get the approximate equation [7], 

A ~ WH 


II. PRELIMINARIES 

Sometimes matrix decomposition has a benefit that it makes 
to reduce dataset and summarization .For example if we have 
an 1000x50 matrix A this matrix has 50000 entities ,if we 
success to factorize the matrix A into a product of two 
matrices W and H such that W is an 1000 x 2 matrix and H is 
2 x 50 matrix , we note that the sum of the entities in the two 
matrices is 2100 > 

Theorem (1) 

Any square n x n matrix A which was written 
as A = BC suc h that , B is an n x k matrix and C j s an 
k x n matrix k<n then A is singular matrix 
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Since k < n we can add n — k zero columns to the matrix 
B to get a new square n x n matrix IV and we can also add 
n-k zero rows or any n-k rows to get another new square 
nxn matrix H 
Now we have 
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Now also we have A = WH 

Also we have det(w) = 0 ,det(//) = 0 ■ 

So we get det(A) = det(iv)xdet(//) = 0x0 = 0 => A is singular. 


Every m x n matrix A can be written as A=BC where B and 
C are two matrices of sizes m x ^ and ^ x n respectively 

where m — ^ 

Proof 

We can form a none singular matrix D=WH where W is of 
type m x ^ and H of type ^ x m 
Now we have an invertible square matrix D of type m x m 
So we can write 

DA = WHA (!) 

Multiplying both sides with D 1 we get 
d'da = D'WHA 

=> A = {d-'w\hA) (2) 

So by taking B = ( D 'W ) , C = ( HA ) theorem (2) will be 
proved. 

Now in equation (2) we have 

a = {d 1 w\ha ) * 

so substituting for the last A in * by A = ( D l w\HA ) , 
then using associative law of matrix multiplication , we get 
A = {D X W)H{D 1 WH)A 
Again we can substitute for the last A to get 


III. METHODOLOGY 

In this paper we want to factorize a matrix into two matrix or 
more using the inverse matrix, so we make use the above 
theorem 

Now suppose we have an m x n matrix A and we want to 

decomoose it into a nroduct of a certain m x m invertible 

matrix D.then we can write 

A = DD 'A => A = d(d 'c)=DC , where (zr‘c)= C 

So we have a matrix A equal to a certain square matrix D 

multiplied by other matrix C , or we can write 

A = (ad )d = ED w jj ere j n this case D is an n x n matrix, 

since multiplication of matrices is associative. 

Now suppose we have an m x n matrix A which is not square 
and we we want to decompose this matrix into a product of a 

certain m x ^ matrix w and other matrix , in this case we can 
form an invertible m x m matrix D by multiplying w by any 
k x m matrix h( i.e D=wh) 

By this method we can factorize any matrix to a product of 
any matrices we want and this can be done by two ways either 

a) we decompose the main matrix into two matrices and then 
we decompose each one of these two matrices using this easy 
method to get four matrices, then again we can decompose 
each one of these four matrices by the same way, and so on 
then by using the association law of matrix multiplication we 
get the result (at every step of decomposition we use different 
invertible matrices) , or 

b) by writing A = AD \ D \ -A, A where D ‘ is invertible 
n x n ma t r j x and us i n g association law to get 

A = (AZ5, Xa a)-(A-iA )A or an y other combination Theorem 
( 2 ) 


A = [d 'W )h(d 1 W )h(d 1 W )HA , and so on 
A = (d'w]h{d'w]h...{d > W )hA = EH EH . . .EH A (3) 

where E = (/Aivjand obviously (/) 'wjH = I m , 
where I m is the identity matrix of rank in 
Since we have I =/ x / ... x / we can deduce that 

m m m m 

I m = (p'w )h(d'W )h..{d'W )h = EHEH...EH (4) 

So we can also analyze the identity matrix of any rank using 
this method. 


The benefit of the easy approach is that ,if a matrix can be 
decompose to certain matrices 

by using this approach we can get the other factors ,to show 
this suppose we have a matrix 
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and we are given two matrices . 
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We want to find a matrix C such that 


T = ABC . 


We can do this as follows: 


AB = 

1 O 

4 2 

2 

5 

1 

£ M 

0 <N 

= 

3 

17 

4 

14 


“—7 

2 

1 

2 





D ' 


13 13 

17 — 3 


_ 26 26 _ 


Next we have 
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T = DD / = Alii) 7 
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Example .1 

Given A _ 


2 3 4 
5 0 1 


1 1 
2 1 


write A= BC where B is 2 x 4 


-1 -2 
1 1 
1 0 
° 5 

call these four matrices standard matrices , or we can use any 
other four matrices satisfying the conditions, in general for 
every matrix A of type mxn we can write A=BC where B 

is m x k matrix and Cis^ xn one such that k - m by taking 
four standard matrices W,H,D, D 1 where W is an 
in x k matrix and H is k x n matrix. 

IV. CONCLUSION. 


D = WH = 


1 -1 

-1 2 


W 


1 1 1 

2 1 0 


,D = 


2 1 
1 1 


matrix and C is ^ x 5 one. 
Solution By taking 



This paper propose very effective and simple method for 
matrix decomposition by using invertible matrices ,so by easy 
method we can factorize any matrix to certain matrices we 
want and other matrices , although there are many methods for 
matrix factorization ,this easy method will be the simplest one 
of them. Unlike the other approach we can decompose any 
matrix to any number of factors (matrices). Easy approach 
has the benefit that for every in X n matrix A we can take 

four standard matrices W,H,D, D 1 where W is an 
mxk matrix, H is kxn to factorize A. 
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Note 

Since matrix factorization is not unique ,we can use more 
simple matrices W, H , D ,D 1 

Also for writing any 2x5 matrix A, as A=BC where B 
is 2 x 4 matrix and C is 4 X 5 one we can use these four 
matrices 
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